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Abstract 



We present a method for the analysis of singularities of Feynman amplitudes 

^^ based on the Speer sector decomposition of the Schwinger parametric integrals 

r~| combined with the Mellin-Barnes transform. The sector decomposition method 

'~~' is described in some details. 

T-H We suggest the idea of applying the method to the analysis of collinear 

K*" singularities in inclusive QCD cross sections in the mass-less limit regularizing 

S-^ the forward amplitudes by an off-shell choice of the initial particle momenta. 

(y~^ It is shown how the suggested strategy works in the well known case of the 

CN one loop corrections to Deep Inelastic Scattering. 

o 

^ 1 Introduction 

. 1^ The success of Feynman's parton model [1] has driven the search for a field the- 

^ ory formulation of the parton hypothesis. The discovery of asymptotic freedom 

[2] of unbroken non-abelian gauge theories has given a clear indication of the 
field content of the Feynman model. Partons are fermions, quarks, and vector 
bosons interacting through exchanges of vector bosons, gluons. 

However, partons, therefore quarks and gluons, should be bound in a hadronic 
initial state. Therefore a high energy partonic process should be characterized 
by two scales, the high energy (momentum transfer) scale and the binding 
energy scale, which should be less than few hundred MeV. 

If one tries to exploit perturbative field theory in order to compute correc- 
tions to the naive free parton model, one finds typically logarithmic factors in 
the ratio of the two scales, which is about one thousand. These logarithmic 
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factors appear as multipliers of the strong structure constant Ug- Thus the va- 
lidity of perturbation theory becomes less obvious in much the same way as the 
idea of parton independence, that is factorization. The only "simplification" 
induced by the presence of the two above mentioned scales is that at least for 
the quarks up and down, whose mass is few MeV, one can consider quarks as 
massless as gluons are. Therefore the field model should be considered in first 
(free) approximation scale invariant and the radiative corrections should violate 
Bjorken scale invariance due to renormalization. This is indeed what happens 
and the corrections to Bjorken scaling [3] give a strong evidence in favour of 
the non-abelian QCD model. 

However the calculations of radiative corrections together with the presence 
of the above mentioned scales, one of which is associated with the process, the 
other appearing in the hadronic wave function, have posed non-trivial technical 
difficulties which have been overcome thanks to the discovery of dimensional 
regularization [4] [5] and of the corresponding minimal renormalization schemes. 
The simplicity of computation is only one of the advantages which have strongly 
favoured the use of the dimensional method. A second, very important, advan- 
tage is that it doesn't break gauge invariance and applies without problems 
to the, however singular, quark and gluon mass-shell amplitudes. This is very 
important in view of parton factorization hypothesis. The dimensional method 
is consistent with the description of the hadron as a gas of free quarks and glu- 
ons whose distribution function has to be computed on the basis of consistency 
conditions of different scale choices. 

The QCD parton amplitudes are however infrared (IR) singular, this hap- 
pens in particular, in the study of deep inelastic scattering (DIS). For this 
reason a huge amount of work has been done aiming at the analysis of the 
singularities of Feynman diagrams as functions of the kinetic invariants, that 
is, of the invariant squared partial sum of the momenta carried by the exter- 
nal vertices. The study of collinear singularities gives information through the 
Altarelli-Parisi [6] approach on the parton distribution function of e.g. a target 
baryon. Infrared singularities appear in a dimensional regularized mass-shell 
massless parton amplitude as poles in the origin of the complex plane of the 
dimensional regularization parameter e = ^^ . Thus in the dimensional scheme 
IR singularities appear in much the same way as UV singularities do. This im- 
plies that one has to avoid IR-UV singularity mixings by first subtracting UV 
divergences. One should however keep firmly in his mind that, while the UV 
singularities are renormalized by a redefinition of few parameters, IR ones are 
associated with a bad choice of the initial and final states of the process. 

An important result concerning the Feynman graph singularities is that in 
the Schwinger parametric representation [7] [8] the singular parts of a diagram 
are confined in particular sectors of the parametric space. This induces a re- 
markable simplification in the singularity analysis but until recently it has not 
been systematically exploited in calculations. (See however [9] [10] [11] and 



references therein). 

Beyond purely technical aspects, an easy tool for the computation of the in- 
frared singularities of the Feynman amplitudes allows a better physical insight 
into Parton Physics. From the point of view of the physical interpretation, 
one should take into account that partons are identified with bound, and hence 
off-shell quarks and gluons. Notice that considering off-shell amplitudes one 
introduces a kind of IR regularization since the infrared singularities are regu- 
larized by a suitable off-shell (Euclidean) choice of external momenta and are 
therefore well separated from the (dimensionally regularized) UV singularities. 
Hence the idea is to regularize the parton amplitudes choosing off-shell ini- 
tial momenta [12] [13]. The singularities of the amplitudes when the parton 
momenta go on-shell are related to the Altarelli-Parisi splitting functions. 

However considering single parton amplitudes with off-shell initial particle 
states presents a further difficulty. Indeed it is fairly well known that off-shell 
charged particle amplitudes are not "gauge invariant". As a matter of fact, 
strictly speaking the single parton contributions are not physically meaningful 
and hence parton factorization, which is the basis of the parton model, should be 
taken "cum grano sails". This difficulty is overcome by the fact that, contrary 
to generic off-shell amplitudes, their mass-shell singular and finite parts are 
gauge invariant and independent. 

We hope that applying the present idea to a physical situation, such as the 
first order corrections to Deep Inelastic Scattering, will help us to make this 
working hypothesis clear. 

In this work we present the off-shell analysis in the light of few, more or 
less, recent progresses in Feynman graph computations. These are essentially 
based on the extended use of Mellin-Barnes transform [14], aiming at the singu- 
larity analysis, and the sector decomposition of parametric Feynman integrals, 
mentioned above. 

A further relevant aspect of our analysis lies in the evaluation of the coef- 
ficients of collinear singularities for Euclidean values of the kinetic invariants. 
These coefficients, which are analytic functions of the kinetic invariants, are 
analytically continued to the physical region where their limiting values au- 
tomatically give well defined distributions, in our case, for the Altarelli-Parisi 
splitting functions. No problem appears on the border of the physical region 
where it is often spoken of soft singularities. We think that the presented 
example clearly exhibits this fact which is however completely general. 

Notice that a generic not renormalized amplitude might present UV sin- 
gularities which might combine with mass singularities. The presence of UV 
and mass singular contributions might appear in the case of multiloop dia- 
grams with UV divergent subdiagrams. However, once these UV divergences 
are minimally subtracted in the Breitenlohner-Maison [15] scheme they should 
not interfere with the mass singularity analysis and, in particular, they do not 
interfere in the one- loop case. 



2 Feynman amplitudes, Schwinger paramet- 
ric form and Speer-Smirnov sectors 

In this section we give a short description of Speer-Smirnov's sectors. We 
start considering a generic connected ^ Feynman ampHtude in d space-time 
dimensions associated with a diagram with / hnes and L loops (and V = I — 
L + 1 vertices). If the amphtude corresponds in the momentum representation 
to a Ld dimensional integral of the product of a numerator N with momentum 
dimension djsf and / scalar propagators, in the Schwinger parametric form the 
same amplitude is [7] [5]: 

j^I+L-Ld V ^ J dM-2a 

^ ) T a=0 b=0 

where Cy is a suitable coefficient depending on the coupling constants and 
combinatorial factors, the sum over r accounts for the different irreducible 
tensors contributing to Ag{p) , @a,b,T is a homogeneous polynomial of degree 

do = d^ — 2a (2) 

in the components of the momenta p entering into the diagram through its 
external vertices. We denote by [X] is the integer part of the positive real 
number X. The parametric integral factors are: 

where 

• the integration domain is identified with the quotient space of the positive 
sector of the /-dimensional Cartesian space and a positive semi-line, that 
is, with the positive sector of the I — 1 - dimensional projective space. 
Therefore, for example, one can write: 



d/i(/3) 



Y[d/3^5{l-J2P^), (4) 



however this is not the choice we shall use in the following, 

Pcil^) is a homogeneous polynomial of degree L in the /3 variables, 

Dg{I3, p) is a quadratic form in the momenta entering the external vertices 
p, whose coefficients are homogeneous polynomials of degree L + 1 in /3, 



^Notice that a diagram is a set of lines and vertices, thus a connected diagram must contain, for 
any pair of its points, enough lines and vertices to form a continuous path joining them. 



• the exponents Xi^r are integer non negative numbers and satisfy the sum 
rule: J2i ^i,T = L{b — a) — a . 

The construction rules of the integrand, in particular, of the Symanzik functions 
Pg(/3) and Dg{I3,p) together with the determination of the exponents Xi^r can 
be found e.g. in [7] and [5] where one also finds the subtraction rules for the 
UV divergences. In [16] one finds the details and proofs of the ensuing analysis. 

A very general IR power counting theorem states that: in a mass-less the- 
ory, such as QCD, in which fields have positive mass dim,ensions and internal 
vertices have dimension four, an amplitude with space-like external momenta 
has no infrared divergences if the external momenta are non- exceptional, i.e. 
no partial sum vanishes. [17] We shall call a theory of this kind, to which the 
present paper is devoted, IR safe mass-less theory. 

Now we come to the sector decomposition of the integration domain and 
to the choice of the parameters /3 in each sector [8] . This is strictly related to 
the structure of the diagram. We shall try to give an idea of the construction 
concentrating on one-particle irreducible (IP-I) diagrams. 

Some basic definitions are in order. We call irreducible a connected diagram 
which cannot be broken into two connected components deleting a line or a 
vertex, we also consider irreducible the trivial case of a diagram consisting in a 
single line. A connected reducible diagram is naturally reduced to a collection of 
parts (pieces [7]) which are irreducible. Given an irreducible diagram G we call 
link a connected sub-diagram Xq of diagram G which contains all the external 
vertices of G and is minimal, in the sense that it does not contain anymore all 
the external vertices if one omits one of its parts. 

The Speer-Smirnov construction is based on the identification of all the 
singularity families of G (s-families) which are maximal sets of G-sub-diagrams 
which are, either links, or irreducible and which do not overlap, that is, either 
they have no lines in common, or are contained into one another. In particular 
the links are always contained into one another. The whole diagram G is the 
first element of any s-family of G. 

It turns out that a s-family is identified assigning a set C of lines of G, that 
we call complete cut, whose deletion breaks G into two tree sub-diagrams each 
of which contains at least one external vertex. The first step in the construction 
is the choice of an order in C. Then the elements of the s-family are identified 
deleting the lines belonging to C in the chosen order^. The first deletion reduces 
G into the union of a link and a set of irreducible sub-diagrams, the other parts 
of what remains of G after the line deletion. Then one continues reducing 
links into sub- links and further irreducible sub-diagrams. When the last link 
is broken into two components, each containing one or more external vertices 
and these components are decomposed into their irreducible parts, one is left 



^In the existing literature instead of cuts one considers the complementary set of lines which are 
called 2-trees. 



with a set of trivial and non-trivial irreducible sub-diagrams. Now one starts 
deleting further lines of the non-trivial elements of the set, breaking all their 
loops. After the last line of C has been deleted, one is left with a set of tree 
sub-diagrams whose parts are trivially irreducible. The s-family is identified 
with the set of links and, both non-trivial, and trivial irreducible sub-diagrams 
that have been generated during this sequence of deletions and reduction into 
parts. 




Figure 1: An example of a s-family in an irreducible 4- loop diagram. The black points 
are external vertices. The solid line sub-graphs denote the links belonging to the s-family 
which also contains the whole diagram and lines 1 and 7. The complete cut in the chosen 
order is the set of lines: (4,6,2,3,5), the reduced cut is (2,3,4,6). A second example of 
s-family corresponds to the ordered cut (4, 5, 3, 2, 6) and reduced cut (2, 4, 5, 6). The reader 
can easily draw the corresponding s-family which now contains 2 instead of 4 links, not 
considering the whole diagram. 

The reader can convince himself that, 

• a s-family F in G contains / elements, if I is the number of lines in G, 

• the number of lines in a complete cut is equal to L + 1 if, as above, L is 
the number of loops, 

• one can establish a one-to-one correspondence l^ between the elements of 
7 G F and the lines of G. This correspondence is obvious for the trivial 
elements of F. In the case of a non-trivial element j , l^ is identified with 
the first line of the ordered complete cut belonging to 7. 

Assigning to every element 7 G F a variable t-y with < t^ < 1 and setting 
to = 1 and /?/ = Yli^^^p i-y , one defines the Speer-Smirnov parametrization of a 
sector of the / — I dimensional integration domain in Eq.(3). This is apparently 
different from that appearing in Eq.(4), indeed for the first line of the cut one 
has /3i = 1 and hence J2i A ^ 1 • Notice that with the new choice of variables 
each integration domain corresponds to a / — 1 dimensional hypercube in the 
t-f space. The reader can verify that the sectors corresponding to different s- 
families are disjoint and the set of sectors covers the whole integration domain. 
On the basis of the above construction of s-families and sectors the para- 
metric integral appearing in Eq.(3) decomposes into the sum over sectors, and 
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hence s-families F: 



lG,{a,h,r) (P) = Y. ^G,{a,b,r) (^) (^) 



with 



r(^) r^\ - ..2Le-rtr Ld f dn^p ULi P^^'^ 



^G 






In the last expression the measure d^Q and the Symanzik functions Pp and 
Dp depend on F, which is the s-family identifying the particular sector. The 
integration measure of the sector amplitude, is 

dix^P = X{d^i= n t\^'^'^dt„ (7) 

1=2 7eF ,77^G 

where the product runs on all the elements of the s-family, except for the 
diagram G itself, and 1(7) is the number of lines contained in the sub-diagram 

7- 

In order to establish the structure of the Symanzik function Pp we have to 
introduce a further kind of sub-sets of lines of G that we call reduced cuts and 
label by C. A reduced cut is the complement of a maximal tree sub-diagram 
of G, it is a set of L lines whose deletion reduces G to a tree diagram. Once 
identified all the reduced cuts, one has: 

PFih) = E n *?""^^ ' (8) 

C I^P' 

where 1(7 n C) is the number of the lines in the intersection of the element 7 
of the s-family and the reduced cut C. 

Let us now consider Dp. On account of the external momentum flow 
through the diagram, given a complete cut C, which does not necessarily co- 
incide with Gp, one defines, up to a sign, the external momentum pc crossing 
the complete cut. Then one has: 

Dp{p,t,) = Y^plY[tl^('rr^^K (9) 

c 7eF 

where 1(7 n G) is the number of the lines in the intersection of the element of 
the s-family 7 and the complete cut G. 

Now we discuss the behavior of the Symanzik functions in a given sector F. 
Since both Pp and Dp are polynomials in the t^ variables and since < t-y < 1 , 
it is possible to identify, given the s-family F, the larger monomial in Pp and 
the larger coefficients of the kinetic invariants in Dp. 

We consider first Pp given in Eq.(8) as a sum over the reduced cuts and we 
single out the dominant contribution. This corresponds to the reduced cut Cp 
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which is obtained from the ordered complete cut Cp identifying the s-family F 
omitting the hne whose deletion breaks the last link. This line belongs to all 
the links contained in F and gives a maximal tree sub-diagram when added to 
the 2-tree diagram made of the lines of G not belonging to Cp] this is not true 
for the ensuing lines in the ordered cut. The monomial corresponding to Cp in 

Pp is given by n76_F^7 ^ • This is larger than the other monomials inside 
the hypercubic sector. Indeed we have: 



Ppih 



7GF C^Cf 'y^^ 



(10) 



with 5^ p(7) = /(7 DC)- I{-i n C^) > and E-yeF %,f(7) > • 

More important for the IR singularities is the corresponding decomposition 
of Dp. In order to analyze the behavior of this function in the sector correspond- 
ing to the s-family F, let us recall, first of all, that one has to study separately 
the coefficients of the various kinetic invariants. Since F corresponds to a given 
complete cut Cp with a chosen order, the corresponding kinetic invariant p"^ 
naturally plays a particular role in our analysis. However its coefficient is in 
general the sum of several monomials, indeed, same momentum pcp crosses 
several complete cuts. If we denote by Cp the set of these cuts, each of which 
gives a contribution equal to n7eF ^7 to the coefficient of the same kinetic 

invariant p^ , one can show that in the sector F the contribution from the cut 
Cp is dominant in the sense that: 



Dp{p,t^ 






CeCT,Cy^CF ■yeF 



fill 






n 4^^'^''"^ 



yeF 



VCF^F{h)+Mp{p,t^] 



where 5c,p{i) = /(t n C) - 1(7 n Cf) > and E^eF ^C,f{i) > • 

Combining our results together we have for the un-subtracted sector ampli- 
tude introduced in Eq.(6) the following expression: 



jiP) („\ _ ..2Le 



r(/ - f - a) dAS!v) 



pI Np{t^) + Mp{p,t^ 



l-^-a ' 



fl2) 



where: 



"'f^G,{a,b,r) 



n 



7eF ,77^G '-7 



^G,la,b,T)'^'' ^ 



dt-. 



(i + E,,c,n,..<^'")^--'' 



(13) 



and 



47aA.)(7) = E ^^ + ^(7) + (^2 " ^ + "^^^^^ "" ^^^ 



-i{jnCF)) + {a-b-^)i{^nCF) . (14) 

Now we see that in an IR safe theory and for Euclidean non-exceptional 
external momenta, that is, all p^ < and d = 4, the amplitude I^ ! ^^ ^n {p) 
diverges, either if / — 2 — o ^ 0, and hence it is a primitively UV divergent, 

or one or more EqIJ^^J'^) < 0, and hence G contains sub-diagrams which 
are UV divergent in the sector F. Indeed, as mentioned above, the amplitude 
cannot be affected with IR divergences in the Euclidean-non-exceptional do- 
main and in this domain the denominator in the integrand in Eq.(12) does not 
vanish. The above divergences must be cured using the Breitenlohner-Maison 
subtraction method. Applying the above formulae to the example presented in 
figure (1), and in particular considering the sectors corresponding to the two 
presented s-families, the reader can easily verify that in the first sector one 
has a contribution which is only quadratically and superficially UV divergent, 
while that from the second sector also presents UV divergences associated with 
the sub-diagram made of lines 1 and 2 and that made of lines 6 and 7. Thus 
the subtraction procedures in the two sectors must be different. We shall not 
discuss this point anymore since in the examples we are going to discuss only 
primitive divergent terms appear which do not contribute to the IR collinear 
singularities we are studying. Thus in particular in our examples the measure 
djl^ (a b t) ^^ integrable in the sector F. 

In order to complete our study we have to change our point of view. Instead 
of looking at a given sector considering the IR singularities developed by the 
contribution of the sector to the amplitude, we must consider one, or more than 
one, kinetic invariants whose vanishing together characterizes the IR singularity. 
We have seen that once the sector is given, this identifies a complete cut and 
hence a kinetic invariant coinciding with the square momentum crossing the 
cut, in general a kinetic invariant corresponds to more than one sector and 
hence several sectors contribute to the IR singularities associated with a set of 
vanishing kinetic invariants. 

Let us consider a set of kinetic invariants with the same (or possibly pro- 
portional) negative value — .^^ and let J-(^ be the set of s-families F^ , and C^ that 
of the corresponding complete cuts Cp. , such that the kinetic invariants |p^ | 

belong to the above set and hence are equal to or smaller than ^^. With these 
definitions the ^^ — )■ O"*" IR singularity appears in the sum over the sectors F^ 



with: 



tCc) („)_„2i. f^„Ft.d r(f - - - a) 



where 

NF,m) = i+ E n4"'"^^'^ (17) 

is a positive polynomial and 

MF,iP,h) = Y. pIU h"^"'^"^ (18) 

where, following our strategy, all the p^ are taken negative. Therefore Mp is 
a negative polynomial which in general vanishes on the boundary points of the 
hyper cube. 

Our purpose is to analyze the behavior of Iq^(^ ^ ^\ when S,'^ — )■ 0+ and this 
can be done employing the Mellin-Barnes transform ^. 

Formally one can write: 



^G.iaM.r^ = 7^ dfl2}:.^, {-MpAp,t,)T^ 



'G,(a,b,T) ~ 2TdJ ^G''("'^'^) ^ ■^« 



2 



X 



Ld , ( S.'^NF^m, 



daT{a)T{I- — -a-a)\- '^ ' \ (19) 

V[-ioo,+ioo\ 2 \^ Mp^{p,t-y)J 

where the path V [—ioo, +ioo] is a continuous line going from —ioo to ioo 
leaving the poles of T{a) on its left-hand side and the poles of r(/ — ^ — a — a) 
on its right-hand side. This path should be closed around the poles of T{a) 
if \C'^Npc{t.y)\ < \Mp({p,t^)\ and around the poles of r(I - ^ - a - ct) if 
\£,'^Np!i{t^)\ > \Mp^{p,t^)\. However, in order to profit of the Mellin-Barnes 
formula, we have to change the order of the integrals considering 

,Fe _ /■ da ,,2.^„^^ .^.. Ld 



,-Fc,d / ,^ / , .^a+^~I ( ^Fi;iP^h)\ /on\ 

df^G,ia,b,r) (-^pMh)) ' I - jv^^(^^) I (20) 

with suitable condition on how the path should be closed. Taking into ac- 
count that Np^it^) is positive and bounded in the hypercube (0 < t^ < 1), if 
—Mp{p,t^) were bounded from below by a positive number, for small enough 
^^ the two formulae in Eq.(19) and in Eq.(20) would coincide with the path 

•^An analogous procedure, aiming at the analyses of power expansions of Feynman amplitudes, 
has been introduced by Pilipp [18]. 
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closed around the poles of r(cr). This is however not always the case, since 
MpApit^) might vanish on the sector boundary. 

A consequence of this vanishing is the fact that the integral 






f211 



instead of being an analytic function of a, might develop a finite number of 
poles. 

It follows that, while in the analytic case the path in Eq.(20) should be 
closed around the poles of r{a), in the second case the path must leave not 
only the poles of T{a), but also those of /(o"), on its left-hand side and those of 
r(/ — 2 — a — a) on its right-hand side and the amplitude must be computed 
closing the path around the poles of T(a) and those of f{(j). 

Let us make this more clear presenting a trivial example: 



Jo C^ + x Jv[-ioo,+ioo] 27ri Jo 



dxx 



a-l 




= -iog(a-f:M^e 

, n! 

n=l 

The pole from the x-integral transforms that of r(o") in cr = from the first to 
the second order and hence produces the term — log(^^) which gives the singular 
part of the original integral. The remaining power series gives log(l + i^^) which 
is the analytic part of the integral around the origin. 



3 Collinear singularities in DIS 

We shall use the lecture notes [19] as a basic reference using the same notation. 
In the framework of the naive parton model the hadronic tensor is computed 
in terms of the parton (quark) distribution function qo{y) and of the partonic 
tensor Wjj^^^{k,q) using the formula: 

W^Ak,q)= ['dyqo{y)Wl,{yk,q), (22) 

Jo 

where y is the fraction of proton momentum taken by the quark. 

In QCD the naive partonic tensor is computed from the tree approximation 
forward quark-virtual-photon Compton scattering amplitude by means of: 
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Thus, according to the optical theorem, the partonic tensor is related to the 
absorptive part of the tree approximation forward Compton amplitude in figure 

(2). 



q+k 



Figure 2: Tree diagram. 

Taking into account the radiative corrections, one has to consider the for- 
ward virtual photon-quark scattering amplitude corrections corresponding to 
the diagrams in figure (3). 

In all these diagrams the initial and final states coincide and consist of a 
quark with momentum k and a virtual photon with momentum q. 

An elastic 2—7-2 scattering amplitude depends on six independent kinetic 
invariants. Four of them are the square momenta (masses) of the four external 
legs, that is of initial and final, possibly virtual, particles. Further invariants 
are the square center of mass energy, i.e. the Mandelstam variable S, and the 
square momentum transfer, i.e. the Mandelstam variable T. 

We denote the independent kinetic invariants by /c^, the square momentum 
of the quark, q^, the square momentum of the virtual photon, S = (q + k)'^ and 
of course T = 0. 

In the case of off-shell quark-virtual photon forward (Compton) amplitude, 
we have negative virtual photon square momentum (g^), positive center of mass 
energy S and null T, while we choose negative virtual quark square momentum 

In principle one should foresee a singularity in the amplitude at T = 0. 
However this is not the case for helicity/angular momentum reasons. A vanish- 
ing angular momentum in the crossed (T) channel Q' + g — )• 7* + 7* is excluded 
by helicity conservation and consequently the amplitude is regular at T = 0. 
The amplitude is not analytic for S* > 0, as expected, indeed we are studying 
the absorptive part of the amplitude. 

However if we complete our off-shell choice considering an unphysical nega- 
tive S together with the negative square momenta of the virtual scattering par- 
ticles, we expect analyticity for infrared power counting [17] reasons. Therefore 
the general idea is the following. 

• We start computing the amplitudes with negative k'^, q^ and 5, where 
they are analytic. 
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Figure 3: One loop corrections. 



• We separate the singular part when k'^ — )■ 0^. The amplitudes are ex- 
pected to diverge as log(— A;^) and in fact they do. We compute the coeffi- 
cient of log(— /c^) in the Mellin-Barnes expansion of each amplitude. The 
Mellin-Barnes expansion gives a complete /c^-expansion of the amplitude, 
however in the present analysis we disregard the contributions which are 
regular in /c^ = 0. 

• The coefficient of log(— /c^) is expected to correspond to an analytic func- 
tion of S, in the region of negative S and q^. By the Mellin-Barnes formula 
and sector decomposition we compute explicitly the analytic coefficient of 
log(— A;^) (as a matter of fact it is the sum of analytic contributions) and 
we analytically continue it in the Re{S) > region, where we find the 
expected branch-cut. We compute the discontinuity, which is directly re- 
lated to the singular part of the trace of the partonic tensor under study. 

We repeat this analysis for each graph, noticing however that, on account of 
the LSZ formula, diagrams d and e have to be divided by 2. 

For each spinorial amplitude we compute the trace of the spinorial matrix 
multiplied by/(f, since we want to sum over the helicities in the (logarithmically 
singular) mass-shell limit. 

In momentum representation and in the Feynman gauge the box amplitude 
is 



A^ul{k,q) = iase^CF 



dp Tr^^P]')^^ 



(24) 



(27r)4 {p-iY{q+pf{k-pY ' 
where cp is the quadratic Casimir of the gauge group. The box contribution 



13 



to the hadronic tensor is: 

WP:i-Hk,q)=Im{A(-;l{k,q)). (25) 

Let's start computing the cohinear divergent box contribution to the trace 
(in Lorentz indices /U and i^) of the partonic tensor. The trace of the hadronic 
tensor gives us the radiative corrections to a hnear combination [19] of structure 
functions. The radiative corrections to an independent Hnear combination are 
given contracting the expression in Eq.(24) with k^k'^. We shall see that this 
contraction does not contain singular contributions in /c^ — )■ 0~. 

Computing the trace we get: 

TrliY^Y<^+:^hu^lp] = IQ {2p ■ kp ■ q + p'^{p ■ k - k ■ q)y (26) 

At this point we pass to the Schwinger parametric form following the procedure 
shown in the former section and disregarding the terms proportional to k"^ in 
the numerator since these give vanishing contribution in the A;^ — t- limit. 
We get 



with K = ase^cpjii^ . Then, integrating over the scale factor t and hence 
passing from the a to the /3-parameters, we get: 



4"'-"-/^(-^ 



2/32 



(/3);d„(p,/3) 



(28) 



with 



Pa{^) = /3l + /32 + /33 + /34 (29) 

We have shown explicitly the infinitesimal imaginary term irj, which ac- 
counts for the time-ordering in the Feynman amplitudes. The amplitude is 
the sum of different terms corresponding to different Speer sectors. We 
focus on the subset of sectors in which the amplitude of the box has 
collinear singularities. These are 6 sectors corresponding to the complete 
cuts crossed by the quark momenta and by the (vanishing) momentum 
transfer. They correspond to the parametrizations: 

a) f3i = s or sf3, (32 = sa, (3^ = sf3 or s, (3^ = 1 ; 

h) /3i = 1 or s/3, (32 = s a, (3^ = sf3 or 1, (3^ = s ; (30) 

c) /3i = 1 or s, /32 = sa, Pz = s or 1, (3^ = s[3 . 
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In this equation each hne corresponds to two sectors, giving the same 
contribution because of the symmetry of Eq.(28) for the exchange of the 
parameters /3i and (3^ . 

Once selected the interesting sectors, we consider again the expres- 
sion for the amphtude given in Eq.(28), which is the sum of two terms. 
The second term, which is proportional to /32, does not contribute to the 
collinear divergence, since /32 vanishes together with Df{p, 13) for k"^ = 0. 

Thus the singularity is only due to the first term in the integrand of 
Eq.(28) and hence we must evaluate the sum of the contributions to 

from the six sectors a), b) and c), that is: 

/iKi) = oKh r s ds da d/3 / 1 

^ "^io {l + s{l + a + ^))^\{l + P){k^+q''sa) + Sa + ir, 

1 1 



(1 + (3 s){k'^ + q"^ a) + Ssa + ir] {I + s){k'^ (3 + q^ a) + S saf3 + i7], 

(32) 

In order to single out the collinear divergence we apply the Mellin-Barnes 
formula, Eq.(20), which gives: 






-2Kk-qf |^(_P)-r(a)r(l-a) (33) 

JT>[-ioD,+ioD] Zm 

sdsdf3da ( {I + (3)-" 



+ 



(1 + s(l + a + (3)Y \ [-«(g2(i + p)s + S) - irf\ 

(l + /3g)-" (1 + g)-" \ 

[-a{q\l + (3 s) + S s) - irf^^-'' [-a{q\l + s) + S s (3) - irf^^''' ) 

Let us perform first the a integral; this has the form: 

'^ da „ ., , 
— a" A{a) , 

OL 

where A[ol) = !/[! + s(\ + /3 + a)] is analytic and non-vanishing in the 
integration domain. It follows that the above integral is equal to A(0)/o" + 
7^((t), where 7^(o") is analytic in the whole complex o"-plane. 

The pole in a = is the singularity considered in the former section 
which must be enclosed in the path encircling the r(o") poles. Thus the 
(T-integrand has a double pole in a = 0, while the pole would have been 
simple in the absence of collinear singularities. 
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Notice that using dimensional regularization as IR-regularization, the 
new pole would be in o" = e. In this case, performing the a-integral one 
would obtain a term — 1/e from the F-pole in a = and {—k'^)~'^ /e from 
the new pole. Then one could "go to the limit" A;^ — )■ 0. Choosing, 
arbitrarily, e real and negative, the second term would vanish in the limit, 
leaving a 1/e singularity, as a memory of the collinear one. 

Following our strategy we compute the cr-integral in Eq.(33), and we 
select the contribution of the double pole since the rest gives analytical 
contribution in k'^ ^ 0. We get: 



4(a) 

J±rp 



r^ 2 /c ■ g , , ,2. /"i sdsd(3 ,^^, 

= K -^ log(-A;^) / 7 , — S;vT^ (34) 

1 1 1 

+ . . . . . - + 



This is an analytic function of the variable S/q^ with a branch-cut in the 
negative real axis. We compute the discontinuity and we set 

4 = 1--' (35) 

q^ X 

where the new variable < x < 1 is the Bjorken variable. 
Setting z = 1 + s(l + a + /3) we find the discontinuity: 



Wt 



p.(«) 



div 



fl fl+2s fjz 

nKlogi-k^) ds - (36) 

JO Jl+s Z"^ 

[5{x z - I) + 5{x z - s) + 5{-z{l -x) + 1 + s)] . 



Notice that —2k ■ q/q^ = 1/x and hence S = q^ + 2k ■ q = g^(l — 1/x). 
Computing the integral factor we get 

^(1 - 2x)e{?,x - 1) + e{2x - 1)^(1 - x) 

+x e{x)e{l - 3a;) + (1 - 2x) e{l - 2x)e{?,x - 1) 
1 — 3t 

+^^e{x)e(i-?>x) , (37) 

where different lines correspond to different sectors and to different phys- 
ical regions. 

Finally the sum of all the contributions gives: 



iyp,(«) 



'T 



log(-^ ) — 7, — (3^ 



div 7f 
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which is the box contribution to the coUinear divergent imaginary part of 
the amphtude we were looking for. 

Before studying the other diagrams, we still have to contract the box 
amplitude Eq.(24) with k^k"^. The numerator in momentum representa- 
tion is proportional to: 

= -32{k-py{k-p + k-q) . 



Trl^Y 






One has 



FFA^fj = -2^^{k-q)'Jdfi{P) 



PaWDa{p,py 



(39) 



(40) 



A'' 



' 2-K 



I 



It is very easy to verify that the amplitude is not divergent in the limit 
/c^ — )■ 0. Indeed, in analogy with the above case we have a /3|-proportional 
term and /3| vanishes with Da{p,/3Y computed for k"^ = 0. 

Considering the contribution from diagram b we have to compute the 
first order in a^ vertex correction which is given by: 

{h}t-hi)l^{i (/3i+/33)+^/33) 
(fc2/3i/33 + g2/3i/?2 + 5/?2/33 + ir/) 

Then we can write the amplitude of diagram b in terms of A'^: 



A^^Ak^q) = -i 



A-K 



S 



17] 



(41) 



In much the same way as for the box amplitude, we first contract this 
tensor with the metric tensor and then with k^k^ . 

Contracting Eq.(41) with the metric tensor, multiplying by two, since 
the diagrams in the figure give the same contribution, and finally disre- 
garding the terms which are proportional to A;^, we have: 



2A 






(/3i + /33 



(5-g2)/33-(5 + g5/32^ 



D,{p,(3) 



2, ,2' 



-log 



A(/3)V 
D,{p,P) 



(42) 



with 



P,(/3) = /3i + /32 + /33 
D,{p,(5) = k^p,p, + q^p,p, + SP2P3 + tV- 
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The sectors contributing to the colhnear singularity are those correspond- 
ing to the parametrizations: 

a) f3i = s, (32 = sa, (33 = 1 ; 
h) (3i = l, (32 = sa, (33 = s . 

In Eq.(42) we have put into evidence two terms. The second one conies 
from the subtracted UV-divergence. It is collinear finite, since it is a 
parametric integral of the logarithm of a function which vanishes on the 
boundary of the integration domain in the k'^ — ?> 0^ limit. The only 
contribution which is collinear divergent in the sectors a) and b), comes 
from the first term, that is: 

24''" = A-|^ / ^(A + A)^|^ . (43) 

S + 17] J Pb{(3) Db{p, (3) 

the rest is finite since it is proportional to (32, which vanishes together 
with Df,{p, (3) in much the same way as for the box. 

Summing the contributions from the two sectors, we get: 

,(6,1) _ 7^2A>g 2n f {l + s)dsda 



2A'^^'' = -K^;^^{S-q 



S + iT]' ^ ' J (l + s(l + a))3 
1 s 



—k"^ + a{—q'^ s — S) — it] —k"^ + a{—q'^ — S s) — irj 



(44) 



We apply the Mellin-Barnes transform and we consider the collinear di- 
vergent part. 
This is given by: 



Z,j\rp 



div 



^^2k ■ q 5, , / ,2n /"^ ds 

K^—^{l-^)\og{-k^) —— 
b +iri q^ Jo (1 + 5, 

(45) 



S + 17] g2 Jq (^i _|_ g^ 

1 s 

s + ^ -ii] 1 + ^ s -it] ^ 



Again this is an analytic function of the variable S/q'^ with a branch-cut 
on the negative real axis and a pole in —irj. 

The novelty of the vertex correction contribution lies in the presence 
of the pole superimposed on the branch cut. Thus computing the ab- 
sorptive part we should take into account the branch-cut discontinuity 
and the contribution proportional to —i7r6{S) coming from the pole. This 
corresponds to a single parton final state which is also present in the 
coUinear-divergent contributions from diagrams d and e. 



However the term i7rS{S) multiplying the vertex correction in diagrams 
b and c is ill-defined, since we see from Eq.(45) that the coefficient of the 
Dirac delta contains the ill-defined integral /q (is/[s(l + s)^] . Furthermore 
the integral factor in Eq. (45) has a branch cut whose discontinuity diverges 
when S vanishes. 

A general remark is here in order. Even if the Fourier transformed 
Feynman amplitude corresponding to a given diagram is the product of 
the contributions from the parts of the diagram, this product structure 
is not suitable for spectral analyses. Indeed for this purpose one has to 
introduce a spectral representation for the whole amplitude and compute, 
e.g. its imaginary part, on the basis of this representation. 

In the light of this comment it is easy to verify that the two above 
mentioned inconsistencies compensate each other. Indeed let us look more 
carefully at Eq.(45). It can be written as follows: 



2A 



div 



K 



logf— A;^ 



X 

1 



ds 



s [l + sf \s+M,-ir^ 



^-^v. 



S/q^ - irj Jo (1 + s) 



1 sds 



\2 S 



1 ds 



1 



s + ^ ~ir] 1 + s-^ -ir] 



(46) 



Now, considering the imaginary part and setting 



S 



i-i, 

X 



(47) 



we have: 



2W^'^'^ 



div 



1 dt 



t{S{x-t)-S{x-l)) 



K 71 log(— A;^) X 

^^ ' h/2l-t 

+ f dt( 6(x - (1 - t)) + 6(x - t)] 

Jl/2 \1 — X J 

-% - 1) (log(2) - i) , 

where we have changed the integration variable t = 1/(1 + s). 

One might wonder if this is a distribution in < x < 1. To check this 
point, we multiply 2W^'^^^'^ \^-^ by a C^^ function 0(a;) and integrate over 

X. 
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We get: 



dx(l){x) 21^^'^^^ 



div 



+ 



Kit log(-A;^ 

1/2 



X 



dx 

1/2 I — X 



(0(x) - 0(1)) 



rfx^0(x)- (log 2-^)0(1) 



Kvr log(— A;^ 



1 , x(t){x) — (f) 

dx 

1 — X 



+ 0(1) 



(4J 



Defining as in the literature [19] the distribution ,-^\ as 



(1 — x)_ 



{l~x)+ 

1 — X 



(49) 



we identify Eq.(48) with: 



2Wr 



T 



as e cp 



lod-k' 



X 



We sum this result with Eq.(38) getting 



+ 6{l-x) 



(50) 



rP,ia) 



{W^'''" + 2W^- 



P,ib)^ 



as e Cp 



div 



log(-A;2 



1 +x^ 



2tt "°' ■" ' \{l-x)^ 
Now one can compute the splitting function P(a;)[19] 



+ 26{l-x) 



/ 1+x^ 
Pix) = cf\ , , +25{l-x) + C5{l-x) 

\(l-x)+ 



(51) 



(52) 



where C is the coefficient of the contribution from the diagrams d and e. 
Notice that diagram f does not contribute to the collinear-divergent part 
of g°''^Wa 13 since it is proportional to the tree diagram term multiplied by 
log(5). 

Considering that, on account of the LSZ formula, diagrams d and e 
have to be divided by 2, one has 



iPXd) 



{W^''"" + w^- 



P,(e)^ 



1 0:5 e'^cp 



\og{-k'^)x6{l-x) . (53) 



div 4 TT 

Hence C = —\ and we get the known value of the splitting function: 

P(a;)=c^(^^|^ + ^5(l-x)j . (54) 
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We have still to contract 2A^^^^ with k^k^. One easily sees that the result 
vanishes in the collinear limit. 

A final remark is here in order concerning the gauge independence of 
our results [13]. This can be proved extending the basic argument for the 
gauge independence of the charged particle mass shell amplitudes in mas- 
sive QED [20]. Indeed in any gauge theory there is a general connection 
between the gauge independence of the "physical" amplitudes and Ward- 
Slavnov- Taylor identities [21]. In the case of amplitudes involving only 
local physical (BRS invariant) operators one finds directly that their par- 
tial derivative with respect to any gauge fixing parameter vanishes. The 
case of amplitudes involving charged fields, which are not BRS invariant, 
is more difficult, indeed one finds that the same derivatives are less singu- 
lar on the mass-shell than the original amplitudes; poles are replaced by 
branch cuts. The idea which can be verified in a general situation is that 
in mass-less theories the same thing happens since any gauge parameter 
derivative of a collinear divergent amplitude is less singular, or even null 
in the mass-shell limit. This is fairly clear in our DIS example, indeed, 
using Slavnov- Taylor identities, one shows that the Feynman-gauge pa- 
rameter derivative of the non-amputated forward scattering amplitude is 
proportional to the box diagram amplitude in which the gluon propaga- 
tor has been replaced by a scalar propagator. The result is apparently 
proportional to the box amplitude (there is a —2 factor), therefore it is 
proportional to log(— /c^). However in order to get the derivative of the for- 
ward scattering amplitude one has to amputate the amplitude multiplying 
it by /c^. This kills the logarithmic collinear divergence proving that the 
gauge derivative of the amputated amplitude vanishes in the mass-shell 
limit. 



4 Conclusions 

We have seen how one can study collinear divergences in DIS cross sec- 
tions by considering the quark mass-shell limit of Feynman amplitudes 
and exploiting the simplifying power of the Speer-Smirnov sector decom- 
position of the parametric integral expressions for the amplitudes together 
with the Mellin-Barnes transform. 

In particular, we have shown that the off-shell regularization is a natu- 
ral tool for the analysis of collinear divergences in inclusive cross sections, 
and that the use of this regularization is strongly simplified by the Speer- 
Smirnov sector decomposition. The recourse to Mellin-Barnes formula is 
the obvious final tool for the computation of the singular parts of the 
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amplitudes. 

The suggestion of initial state off-shell regularization of inclusive cross 
sections as a reasonable "physical option", as an alternative to the widely 
adopted dimensional scheme, is one of the results of the present paper. 
However the principal goal of this paper was to identify a general method 
of analysis of mass singularities in the Feynman amplitudes of massless 
field theories, which could be automatically applied to multi-loop and 
multi-external-vertex diagrams. For this reason we have described with 
great care the role of different kinds of cuts in the construction of the 
parametric form of a Feynman integral. These cuts are sets of lines in 
Feynman diagrams whose deletion either breaks the original diagram in 
two parts, or transforms it into a tree diagram, or else does both things 
together. We have also analyzed the role of the cuts in the identification 
of the Speer-Smirnov s-families and sectors. The idea is that the iden- 
tification of sectors contributing to the IR-singular parts associated with 
the vanishing of kinetic invariants is a central step in the study of these 
singularities, and that the method could be extended to multi-loop and 
multi-external-vertex diagrams developing suitable software tools. As a 
matter of fact, in the last year an example of such tools (Fiesta 2 [11]) 
was developed by Smirnov and Tentyukov. However, to our knowledge, 
the method has not yet been used in the explicit analysis of collinear 
divergences. 

A systematic use of the techniques presented in this paper may in 
principle prove useful in different contexts of phenomenological interest. 
High-order perturbative calculations in QCD, for example, may take ad- 
vantage of a systematic separation of singular contributions. From a dif- 
ferent point of view, an efficient way of identifying IR divergent terms to 
all orders in perturbation theory can be useful when one is faced with the 
problem of resumming the whole perturbative series in special kinematic 
regimes, where powers of large logarithms of ratios of invariants spoil the 
reliability of fixed-order calculations. 

A point that we have not discussed with sufficient completeness is the 
effect of the UV subtractions parts on the IR analysis. A further study 
on this point is crucial if one tries to work with high-order diagrams. 
One should compare Speer-Smirnov s-families with Hepp-Zimmermann 
forests and show that the UV subtraction should be limited to the sec- 
tors whose corresponding s-families contain one-particle-irreducible UV 
divergent components. 

In conclusion, being clear that Speer's sector decomposition is useful 
in the singularity analysis, while it does not seem particularly suited for 
regular diagram calculations for which other sector decompositions have 
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been developed together with a suitable software [18] [9] [10] , we believe 
that our work has been convincing enough concerning the advantages of 
the Speer-Smirnov sector decomposition combined with the Mellin-Barnes 
formula in the IR- singularity analysis, also in view of the study of multi- 
loop multi-leg amplitudes. This point is also stressed by A.V. Smirnov, 
V.A. Smirnov and M. Tentyukov in some recent works [22] [11]. One 
should however notice that, when the number of legs increases together 
with that of the corresponding kinetic invariants, the analytic continua- 
tion foreseen in our strategy might become prohibitive. It is clear that 
the discussed example and the Drell-Yan processes, in which the analytic 
continuation concerns a single variable, the total energy, correspond to 
particularly lucky situations. 
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